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Abstract 

The issue of the local visibility of the shell-focussing singularity in 
marginally bound spherical dust collapse is considered from the point 
of view of the existence of future-directed null geodesies with angular 
momentum which emanate from the singularity. The initial data (i.e. 
the initial density profile) at the onset of collapse is taken to be of 
class C 3 . Simple necessary and sufficient conditions for the existence 
of a naked singularity are derived in terms of the data. It is shown 
that there exist future-directed non-radial null geodesies emanating 
from the singularity if and only if there exist future-directed radial 
null geodesies emanating from the singularity. This result can be 
interpreted as indicating the robustness of previous results on radial 
geodesies, with respect to the presence of angular momentum. 
Pacs:04.20.Dw, 04.20.Ex 



1 Introduction 

It has long been known that naked singularities may arise in the gravitational 
collapse of inhomogeneous dust spheres; see [Q , where the existence of naked 
so-called shell-crossing singularities was first demonstrated. The physical 
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significance of these singularities has been questioned, primarily on the basis 
that they are gravitationally weak according to the definition of Tipler 0, 
[|. However, it remains to be shown how to obtain a unique evolution to 
the future of a shell-crossing singularity. A more disturbing type of naked 
singularity was discovered and christened the shell-focussing singularity in 
Q. This singularity was further studied in || and 0, where the first results 
regarding the role of regular initial data in determining the causal nature 
of the singularity were derived. This question has been studied extensively 
since then (see || and references therein). The main conclusion of these 
studies is that for a given initial mass distribution, the remaining initial 
data (namely the velocity distribution of the dust cloud) may be chosen so 
that the collapse results in either a black hole or a naked singularity |7J. 
Furthermore, there are open subsets (of the appropriate initial data space) 
of choices for the velocity distribution which lead to these conclusions. In 
this sense, both censored and naked singularities are stable [§, ||. Indications 
of the instability of the Cauchy horizon associated with the naked singularity 



have arisen in perturbative studies [10 



The aim of this study is to address the stability of the naked singularity 
from another point of view. All results to date on the visibility of the singular- 
ity relate to the existence of future-pointing radial causal geodesies emanating 
from the singularity. (In all but one case [Til], where radial time-like geodesies 



are treated, the results relate solely to radial null geodesies.) Our goal is to 
investigate the occurrence or otherwise of non-radial null geodesies emanating 
from the singularity. This is of importance as non-radial geodesies provide 
a better model of physically realistic trajectories than do radial geodesies. 
On the mathematical side, radial null geodesies passing through a point p 
of space-time constitute a set of measure zero in the set of null geodesies 
through p and so must be considered to be extremely specialized. Further- 
more, investigating non-radial geodesies gives an indication of how angular 
momentum (at the level of the metric) might influence the system. 

The structure of this paper is as follows. In Section 2, we discuss the 
field equations and occurrence of singularities and give a clear definition of 
the initial data which are being studied. This is essential for a rigorous 
study of the possible outcomes. The initial data consist of just one function 
/j,(r) of a single variable on an interval [0,6], where b is the initial radius 
of the dust cloud. We take /i G C 3 [0,6]; this is more general than the 
C°° or analytic data often studied. In Section 3, we derive necessary and 
sufficient conditions, in terms of the initial data, for the existence of radial 
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null geodesies emanating from the singularity. Our results agree with those of 
previous studies || |7|, [L2| , but have the two main advantages of (i) overcoming 
some potential mathematical deficiencies of previous approaches, for example 
in the case of the proof of necessary conditions for the occurrence of naked 
singularities; (ii) giving a complete decomposition of the C 3 initial data space 
into regions leading to naked and censored singularities. Also, the analysis 
used here is quite different from the approach of previous authors and sheds 
some new light on the nature of the singularity. We study non-radial null 
geodesies in Section 4. Our main conclusion is this: there exist non-radial 
null geodesies emanating from the singularity if and only if there exist radial 
null geodesies emanating from the singularity. Some concluding comments 
are given in Section 5. All our considerations are restricted to the local 
visibility of the singularity. A bullet • indicates the end of a proof. We use 

87TG = c= 1. 



2 Field equations and singularities 

We study spherical inhomogeneous dust collapse for the marginally bound 



case. The line element is [0, [14], [T5[ 



ds 2 = -dt 2 + {R r dr) 2 + R 2 {r, t)dQ 2 , (1) 

where dVL 2 is the line element for the unit 2-sphere. We use subscripts to 
denote partial derivatives. The field equations yield 



R, = -ff , (2) 

p = wt r - (3) 

the latter equation defining the density p(r, t) of the dust; m(r) is arbitrary^ . 
(m will be referred to as the mass, although in fact m equals twice the Misner- 
Sharp mass. In particular, this locates the apparent horizon at R = to.) 
Solving (g) gives 

i? 3 (r,t) = %(t c (r)-t) 2 , (4) 



1 m(r) will, in fact, be taken to satisfy certain differentiability conditions and inequalities 
defined ahead in the text. 
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where t c (r) is another arbitrary function which, in the process of collapse, 
corresponds to the time of arrival of each shell r to the singularity. 



Assumption 1 The collapse proceeds from a regular initial state; i.e. at 
time t = 0, (i) there are no singularities (all curvature invariants are finite) 
and (ii) there are no trapped surfaces: 

m(r) < R(r,0). (5) 

We have the freedom of a coordinate rescaling in (|l|) which we can use to set 
R(r, 0) = r. This gives 

tc(r) = \J^. (6) 

Notice that this leaves just one arbitrary function, m(r). There is a curvature 
singularity called the shell focussing singularity along t = t c (r), so the ranges 
of the coordinates r, t are < r < oo and < t < t c (r). Assuming a 
dust sphere of finite radius, we can restrict the range of r to [0, b] for some 
b > 0, and match it to a Schwarzchild exterior. During the collapse of the 
dust sphere there can be also a curvature singularity given by R r = along 
t = t.Jr), where 



, . wrm , . 

t sc r = 2^ . 7 

m ,r 

It is known that this so-called shell-crossing singularity is gravitationally 
weak fH, H, though what this means in terms of continuability of the geometry 
is not yet known. Since we are primarily interested in the shell focussing 
singularity we impose: 

Assumption 2 Along each world-line r =constant, the shell- crossing singu- 
larity does not precede the shell-focussing singularity. That is, t c {r) < t sc (r) 
for all r > 0. 

This is equivalent to taking R r > for all r > 0, and yields 

j ii 

m r <3— , r > 0. (8) 
r 

We note that a no-shell crossing condition should imply t' c (r) > and, in 
this case, this is in fact equivalent to (|). 
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In order to obtain rigorous results relating to the initial data, we need to 
be very clear about what those initial data are. This can be phrased in terms 
of the class of functions to which m(r) may belong. The freely specifiable 
initial data consists of just one function: the initial density 

p(r):=p(r,0) = ^f. (9) 

Through the next assumptions we ensure the positivity of the matter-density 
p during the collapse. We impose a certain differentiability level on p which 
is higher than one would normally like to assume (C 1 ), but allows for a 
comprehensive discussion of censorship. 

Assumption 3 1. The mass function m(r) is strictly increasing on [0, 6]. 

2. The density function p{r) is such that 

peC 3 [0,b], tio:=n(0)^0. (10) 

We define three functions m i; i — 1, 2, 3 as follows. The definitions are meant 
to apply for whatever values of r the quantities exist. 

m(r) =: r 3 (— + m 1 ), mi(0) = 0; (11) 
3 

«V(r) =: r 2 (p + m 2 ), (12) 

m tr r(r) =: r(2/i + m 3 ), m 3 (0) = 0. (13) 

Notice that m 2 (r) = p(r) — p , and so m 2 € C 3 [0, b] and m 2 (0) = 0. 

We now turn to the question of how the no-shell crossing condition re- 
stricts the initial data m(r). In order to do that we prove a number of results 
which constrain both m(r) and its derivatives. 

Proposition 1 The no-shell crossing condition (W is equivalent to 

m 2 (r) < 3mi(r),r e (0,6], (14) 

which is equivalent to 

m[(r) < 0,r <E (0,6]. (15) 
Proof: This is immediate from the definitions above* 
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mi,m 3 


G C 3 [0,6] and 






mi(0) 


= ^(0), 


m 3 (0) = 


3m' 2 (0); 


m?(0) 


= ^(0), 


m'M = 


4m' 2 '(0); 


m;"(0) 


= ^<(0), 


<(0) = 


= 5<(0) 



Corollary 1 TTie no-shell crossing condition yields m 2 < 3m! < on (0,6]. 

It is important to study in more detail the properties of m(r) at r = 
since they will play an important role in deciding whether a null geodesic can 
escape from the singularity. We will now demonstrate some conditions on 
the values of the derivatives of m(r) at r = 0. The proof of the next result 
makes use of elementary calculus to describe some important relationships 
between the derivatives of the functions raj(r), i = 1,2,3 which will be of 
use in the following section. 



(16) 
(17) 
(18) 

Proof: Clearly, it is only differentiability at the origin which needs to be 
checked. By the definitions above, we can write 

1 f r 

m i( r ) — / x 2 rri2(x) dx, (19) 
V Jo 

1 f r 

1^2(1") = — J xm 3 (x)dx. (20) 
Using Taylor's theorem to first order, we can write 

rri2(x) = xm' 2 (x), (21) 

where the mapping x 1— > x(x) : [0,6] — * [0,x] is continuous. Then 

1 f r 

m i( r ) = ~ x 3 m' 2 (x)dx 
r 6 Jo 

r 
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where in the second line we used the generalised mean value theorem for 
integrals. The mapping r i— > f(r) : [0,6] — > [0, r] is continuous. Then 

mi(r) — mi(0) , mi(r) 
hm = hm 

r— >o+ r r-»0+ r 

= lim -m'Jf) 

r ^o+ 4 

= lirn -m^(0), 

i — >0+ 4 

which says that m! x (0) exists and equals the quantity asserted. The results on 
the higher derivatives of m\ are obtained in a similar manner. Differentiating 
(|20D and rearranging, we can write 

-fl = m' 2 (r) + 2^1. (22) 

then taking the limit and using l'Hopital's rule for the second term on the 
right hand side shows that m' 3 (0) exists and equals the quantity asserted. 
Results for the higher derivatives of are obtained by analysing the deriva- 
tives of (|22p • 

As a corollary of the above results one can write useful conditions on the 
derivatives of the functions rri2 such that the no-shell crossing condition is 
satisfied. 



Corollary 2 The no-shell crossing condition yields: 

1. m' 2 (0) < 0. 

2. Ifm' 2 (p) = 0, then m' 2 '(0) < 0. 

3. J/m 2 (0) = m 2 '(0) = 0, then m' 2 "(0) < 0. 



To summarise the main points of this section, we have proved a number of 
results which establish conditions on the mass function such that the process 
of collapse starts from a regular initial state and the shell- crossing of the 
dust spheres is avoided thus ensuring that the collapse will end in a (central) 
shell-focussing singularity. In what follows we will consider the question of 
whether this singularity is naked or covered. 
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3 Radial null geodesies 



We recall that a singularity is called locally naked if in its neighbourhood there 
are future-directed causal geodesies which emanate from the singularity. In 
this section we study the case of radial null geodesies^. Using (P, (|) and 
© one can write the equation governing the outgoing radial null geodesies 
as 

dt „ 1 2 / fix 3 V 3 /„ / \ 



dr 

We know that if a radial null geodesic (RNG) emanates from the shell fo- 
cussing singularity, then it must immediately move into the regular (un- 
trapped) region. That is, it must precede the apparent horizon, which is 
given by R = m, i.e. by 

* = t H (r) = j^r 1/2 - \m > 0. (24) 

Furthermore, the only 'point' of the shell focussing singularity which may be 
visible is r = 0, whereat t c (0) = f (^jf) - ^ 2 =: *o- So we are seeking RNG's 
that extend back to (r = 0, t = to)- Henceforth, we shall refer to this point 
as 'the singularity' which will be called naked if it is, at least, locally naked. 
We shall now derive sufficient conditions on m(r) such that the singularity 
is naked. The main idea of the proof is that in order to find geodesies 7 
escaping from the singularity one must ensure that they lie below a curve 
given by t = tzj(r) (i.e. they precede the time of the apparent horizon) and 
above another curve t = t*(r) which ensures R <rt > as well as the no-shell 
crossing condition. This idea is represented schematically in Figure [1|. 

Proposition 3 If m' 2 (0) 7^ ; or if m' 2 (0) = and m'^O) 7^ 0, then the 
singularity is naked. More precisely, defining 

U(r) = - ^™ > 0, (25) 

r 6 m tr 

then there exists c > such that every outgoing radial null geodesic which 
passes through the region 

Snaked ■= {(r,t) i < t < a, t*(r) < t < t H (r)} 



2 In fact the results of this section apply to arbitrary outgoing radial null curves; an 
affine parameter is not introduced at any stage. 
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originates at the singularity. 



Proof: First, we note that tjj > t* if and only if 

1 2 Q/Z^o \s/2/ \ 
mi > -m 2 + -r { — + mi) 1 (p + m 2 ), 

and so subject to the hypotheses of the Proposition and the conditions on the 
mass, there exists a sufficiently small c so that the region Q na ked is non-empty. 
Note also that t#(0) = £*(0) = to- It is easily established that t#(r),t*(r) 
are differentiable on [0, b], and we find 

Ur) = -±m> 2 (0)(^r 3/2 - ^(OXf )' 3/2 r + 0(r 2 ). (26) 

Now let p be any point of Q, na ked- Then there exists a unique solution of ( |23| ) 
through p, which for r > is a differentiable curve t ra d{r). Decreasing r along 
t ra d, the RNG must either cross t — tn, t — t* or run into the singularity. 
Using (p3|), we find that 

Kad\t=t H = 1 'h + m ,r > t' H for all r G [0, 6], (27) 

with equality holding only at the origin. Thus in the direction of increasing 
r, t ra d crosses from below, and so extending back to r = 0, we see that 
trad must stay below We can also calculate that 

t' rad \ t=tt = 2- 2 / 3 (^ + miY\pi + m 2 y/ 3 (3 mi - m 2 ) 2 / 3 . (28) 

It is then a straightforward matter to check that in each of the two cases 
referred to in the hypotheses, there exists al > such that, for all r e (0, d), 

Kad\t=t t < K( r )- 

Thus in the direction of increasing r, t ra d crosses £* from above, and so 
extending back to r = 0, we see that t ra d must stay above £*. 

Hence we see that the geodesic t ra d must extend back to the singularity 
(r = 0,t = t ) • 

We note that one has t > £*(r) if and only if R tTt > 0. This will be of 
importance when we turn to the study of non-radial null geodesies. We also 
note that the condition £* > gives 

m , , 

m r >—. (29) 
r 
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In the last proposition we have shown sufficient conditions for nakedness. We 
will now show necessary conditions on m(r) such that the singularity is cen- 
sored and which help us to isolate in a final case the remaining possible naked 
solutions. As in the previous proposition the proof follows by comparing the 
slopes of the curves depicted in Figure [I]. 




Figure 1: Schematic diagram of the singularity structure in the radial 
2— space. In the region Q = {(r,t) : r < a,t* < t < tn} no outgoing 
radial null geodesic 7 can cross either tu or {t q in the final case), as r 
decreases, and so must extend back to the singularity (r = 0, t — t Q ). If tn 
is decreasing in a neighbourhood of the origin, then the singularity must be 
censored. 

Proposition 4 lfm' 2 (0) = m'^O) = and (^"(O)! is sufficiently small, then 
the singularity is censored. 

Proof: We know that dt ra d/dr > for r > and that a RNG cannot emanate 
into the trapped region. Thus a necessary condition for the existence of 
a naked singularity is that tn{r) is increasing, i.e. t' H (r) > 0, on some 
neighbourhood (0, a) of the origin. This is equivalent to tuij) > t*(r) on 
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(0, a) , which in turn is equivalent to 

u(r) := mi - ]-m 2 - \r z (^- + mi) 3/2 (/i + m 2 ) > 0, re (0, a). 

Subject to the current hypotheses, we find that 

lim ^ = -2(^) 5 / 2 - i<, 
r^o+ r 3 v 3 ; 36 2 ' 

and so the conclusion follows • 

The following corollary is immediate: 

Corollary 3 If m' 2 (0) = rn^O) = m 2 (0) = 0, then the singularity is cen- 
sored. 

This is a very useful result as it shows that there is only one case left to 
consider, that where m' 2 (0) = m'^O) = but wl^ifS) ^ 0. Thus the structure 
of the mass function is completely decided for the consideration of this case, 
which we call the final case. 

3.1 The final case 

We assume for the remainder of this section that m' 2 (0) = m^O) = and 
m 2"(0) < 0. The sign of the last term here is required by the no-shell crossing 
condition. It is convenient to define the number (3 > by 

<(0) = -l2^f/\ 

Then from Proposition we have, for r — ► + , 

mi(r) ~ -^(y) 5/ V, m 2 (r) ~ -2(3{^f'r\ (30) 

Furthermore, from the proof of Proposition 4, we see that we can restrict 
our attention to (3 > 6; if (5 < 6, then the singularity is censored. We will 
now show that with our approach one can also easily recover the known 
result |F| [L6| that, in the final case, sufficient conditions for the existence of 
naked singularities are obtained through the solutions of an algebraic quartic 
equation. Again the proof follows by calculating the slopes of the curves in 
Figure |l| 
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Proposition 5 Subject to the conditions fffl]) which define the final case, 
sufficient conditions for the singularity to be naked are that (3 > 6 and that 
there exists at least two positive roots q G (0, oo) of the equation 

2g(^) 1/3 = 3/?-g. (31) 

Proof: Let q > and define t q (r) = t + f(^f-)?" 3 - Notice that t q (r) intersects 
the singularity at r = 0. Then the region {(t,r) : < r < a,t q (r) < t < 
tnir)} is non-empty for some a > if 

(3 > q + 6 > 6. 

We calculate that the slope of a radial null geodesic (RNG) crossing t q is 

R r \t=t v (r) = (y)(— g— ) ' (— ) r +°( r )> 



while the slope of i„ is 



,/ / \ Q I t^O \ 2 



Thus the former slope is less than the latter, and so the RNG cannot move 
below t q as r decreases (hence giving a naked singularity), provided 

( ff-g )-i/3( 3 ff-g ) < £ 
6 6 3 

The result follows • 

The converse of this result is also true, but is considerably more difficult to 
prove. We recall that in past works f| 0, [12| either the emphasis was given on 
obtaining only sufficient conditions for the occurrence of a naked singularity, 
or the tangent to a null geodesic emanating from the singularity was assumed 
to have a finite, positive limit in a particular co-ordinate system. We note 
that, in principle, the non-existence of such limit might not be a sufficient 
condition for the non-existence of a geodesic emanating from the singularity. 
We give an example of this in the last section. We show next a sketch of 
the rigorous proof that the sufficient conditions given by Proposition [5] for 
the existence of a naked singularity are in fact also necessary. The details of 
the proof are omitted. Unfortunately, the usual fixed point theorems do not 
apply, and so an ad hoc approach is required. 
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Proposition 6 Subject to the conditions ffiJQ) which define the final case, a 
necessary condition for the singularity to be naked is that there exists a root 
q G [0, oo) of the equation 

2g(^) 1/3 = 3/?-g. 
o 

Sketch of proof: In this final case, the governing equation (|23|) for the 

geodesies may be written as 

d f = 4 )" 1/3 (1 + A(«))(l + / 2 («))- 1/3 (l - \ V ~ + fs(u)), (32) 

where y = 6/3 -1 (^) _:L (i — to), u = r 3 and fi(u),i = 1,2,3 are C° functions of 
m which vanish in the limit u —>■ + . We assume that there exists a solution 
yi(u) of (j3^); i.e. G C 1 (0, Uj] for some Ui, yi(0) = and y\ satisfies (j32|) 
on (0, Ui]. Furthermore, the solution must satisfy, for u G (0, v^], 

< yi(w) < y H (u) =au + / 4 (tt), 

where a = (/3 — 6)//3, y# corresponds to the apparent horizon, and this 
defines the continuous o(it) function f&(u). 

A key point of the proof, which is used repeatedly, is that if the limit 
lim n ^ + u ^Ui exists, then using l'Hopital's rule and (0), this limit must 
satisfy a certain quartic equation which strongly restricts the allowed values 
of this limit. In fact we do not need to prove that this limit exists as the 
quartic equation arises elsewhere in the analysis. Furthermore, the existence 
of this limit is not guaranteed by the proof below, and the limit possibly may 
not exist. 

The proof proceeds by iteratively improving the bounds < yi/u < Vh/u 
by using existing bounds in fl32|) and integrating. Unfortunately this is not 
a straightforward process, as the bounding intervals do not automatically 
contract at each step. However, by studying carefully the subregions of a 
current bounding interval where contraction does not occur, we can show 
that the resulting 'expansion' leads to a contradiction such as yi/u > a. The 
first step is to show that there exists uq > such that 

k < — <l 
u 
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on (0, uq], where k , l are respectively the smaller and greater positive roots 
of go(x) = x 4 — x 3 + Aq, where 

3 _ 6 /3 + 3 3 

Such roots must exist in order that yi exists. Next, we define the sequences 
{A n }^L , {/c n }^o an d {^n}^=o as follows. With the n = terms as above, for 
n > we take 

v, = §(f)"°a-|), 

and define k n ,l n to be the smaller and larger positive roots of g n +i( x ) = 
x 4 — x 3 + A 3 ,!- As for the first step, we can force the contraction u~ x y\ G 
(k n ,l n ). Clearly, k n+ i,l n+ i are not defined for all n for particular values 
of (3. When this happens, u~ l yi leaps out of the current interval (k n , l n ) 
giving a contradiction and so a solution cannot exist. So we assume that (3 
is such that the three sequences are defined. Then it is readily shown that 
{A n }^0' {^n}^=o are strictly increasing, while {l n }^ =0 is strictly decreasing. 
If {A n }^ is unbounded, we obtain a contradiction. Hence this sequence, 
and consequently {/c n }^0' {^«}^o mus t converge. It is easily demonstrated 
that the limit I := lim^oo l n must be positive and must satisfy 

^(i-0 1/3 = |(|) 1/3 (i-^)- 

Taking q — (31 proves the result* 

We note that the proof avoids using the assumption that lim n _ +0 + u^yi 
exists, nor does it conclude that this is the case (i.e. we do not have 
lim^oo k n ^ lim^ooZn). The necessary and sufficient conditions for the 
occurrence of naked singularities in the final case can be summarised in the 
next result which is obtained by carefully studying all the possible roots of 
©: 

Corollary 4 In the final case, excluding (3 = (3* := |(26 + 15\/3) ~ 77.97, 
the singularity is naked if and only if (3 > (3*. 

Proof: We will see below that for (3 = the equation (|3T| ) has exactly one 
positive root. This does not guarantee that the sufficient condition of Propo- 
sition 5 is verified, and so we exclude this case from our considerations. Then 
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as we have seen, the necessary and sufficient conditions for the singularity to 
be naked are that (3 > 6 and that there exist at least two positive roots of 



fx(l) = l 4 ~l 3 



+ A(l--) 3 = 0, 



where A = 81/(4/3). The values of A for which fx has roots are determined as 
follows. First, we look for double roots of fx, i.e. values of / for which f\(l) = 
/a(0 = 0- These occur for two values of A only; A = Ai = 27(26 — 15y / 3)/2 
with double root l x = 3(2 - \/3) and A = A 2 = 27(26 + 15^/2 with double 
root Z 2 = 3(2 + y/3). Notice that since l\ < 3, fx{h) < for all A < Ai, 
and since l 2 > 3, fx(h) < for all A > A 2 . Thus fx has positive roots for 
A G (0,Ai) U (A 2 ,oo). 

Suppose next that fx has two positive roots for A G (Ai, A 2 ) Then f' x has 
a positive root U which can be shown to lie in (/i,/ 2 ). But the conditions 
A G (Ai,A 2 ), f'{U) — 0, G (hjh) yield fx(l*) > 0, which cannot be true. 
Hence fx cannot have two positive roots for AG (Ai, A 2 ). 

Thus fx has (at least) two positive roots iff A G (0, Ai) U (A 2 , 00). In terms 
of 0, this is iff (3 G (0, 3(26 - 15^/2) U (3(26 + 15-^3) /2, 00). The condition 
(3 > 6 restricts us to the latter interval • 

For convenience, we summarise the results of this section as follows, phras- 
ing our results in terms of the behaviour of the initial density function n{r). 

Theorem 1 Let b > and /i G C 3 [0, 6]. Let 



bound collapse of the dust sphere with initial radius b and initial density pro- 
file n(r) results in a naked singularity if and only if one of the following 
conditions is satisfied. 




defined on [0, b] satisfy the no-shell crossing condition (Qj. Then the marginally 



1. //(0) < 0. 



2. n'{0) = and //(0) < 0. 



3. n'(0) = /i"(0) = and 



< -- 




(MO)) 5 ' 2 
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We note that the results of this section are in complete agreement with 
those of previous studies |12|, [16], |17| . Our approach, which will be also used 
in next section, is substantially different and so gives independent verification 
of these past results on radial null geodesies. Furthermore, our treatment of 
the final case is more general than it was previously done since we establish 
rigorously that, from the RNG's analysis, the sufficient conditions for the 
occurrence of naked singularities are also necessary. We will comment on 
this again in the final section. 

Having determined fully and rigorously which data give rise to singular- 
ities which are radially visible, we may now consider the case of non-radial 
null geodesies. 



4 Non-radial null geodesies 

In this section, we determine the circumstances under which a non-radial null 
geodesic may emanate from the singularity. This is potentially important 
since, physically, one can expect null geodesies to have some angular mo- 
mentum. Furthermore, by introducing angular momentum in the geodesies 
equations we will also be doing a stability study of the naked singularity 
solutions that appear in the radial case. 

The governing equations for the non-radial null geodesies (NRNG's) ob- 
tained from the Euler-Lagrange equations are 

R rt 'o L 2 , R + R r t . . . 

t+~ffl 1 +!T2^--f) = °» ( 33 ) 

it T it it it r 

L 2 

R r r + R rr r 2 + 2R rt ri - — = 0, (34) 

R 6 

L 2 

-i 2 + (R >r r) 2 + — = 0. (35) 

The over-dot represents differentiation with respect to an affine parameter s 
along the geodesies, and the constant L is the conserved angular momentum. 
We are looking for is the existence of a solution of (|33| - p5D , which satisfies this 



Existence condition: 

There exists e > such that t is a non-negative, integrable function of the 
affine parameter s and r is an integrable function of s for s G [0, e) and such 
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that 



lim r(s) = 0, lim t(s) = tn. 



Non-negativity of t implies that the geodesic is future-pointing. We im- 
mediately see from (p5|) that t must be infinite in the limit s — * + , and hence 
we must have t —>■ — oo as s —>■ + . In fact this applies to any NRNG which 
extends back to the center R — (r = 0). These conditions will be used 
in the proofs that follow. The main idea of the proofs given in this section 
is essentially the same as in the case of RNG where in order to have naked 
singularity solutions, one needed to find geodesies 7 which lie in the region 
between t H and (as shown in Figure [I]). The existence of the lower bound 
t* obtained in the next proposition is a simple consequence of the geodesic 
equations 



Proposition 7 A non-radial null geodesic emanating from the singularity 
must lie in the region t*(r) < t < tj^(r) for all r G (0, a) and some a > 0. 

Proof: The second inequality is obvious. For the first, note that using the 
geodesic and field equations we can write 



and so 



t>-^f{R r rf. (37) 

it r 



Thus for a solution which extends back to the singularity, we must have 
R t rt > 0, which gives t > t*(r) (recall that this is the defining property of 
t*(r)) • 

We can now consider the following problem: Suppose that for a given 
mass function m(r) there are no RNG's escaping from the singularity. Can 
we find NRNG's which escape? Before giving the answer below in Proposition 
[| we need to prove the following result: 

Lemma 1 r ^ in a neighbourhood of s = + along an NRNG which 
extends back to the centre r = 0. 

Proof: Suppose that f = at some point s = s on such a geodesic. From 



(36) this gives i positive which cannot happen sufficiently close to r = 
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Note that this lemma applies to the case where the geodesic extends back 
to the singularity. We can now give an answer to the previous question 



by simply using the geodesies equations ( |33| - p5|) and the results of the last 
section on RNGs: 

Proposition 8 If the singularity is radially censored, then it is censored. 
That is, if there are no radial null geodesies emanating from the singularity, 
then there are no null geodesies emanating from the singularity. 

Proof: Assume that the singularity is radially censored and suppose that 
there exists an NRNG which extends back to the centre. By Lemma 1, we 
have f 7^ and since generally _R r 7^ 0, we can rewrite as 



t ,0 , L 



l + (-^-) 2 > 1. 



R r f RR T f ' 

Thus along an NRNG, we have 

— > R r . 

ar 

The right hand side here gives the slope of an outgoing radial null geodesic, 
the left hand side is the slope of the NRNG's. Now let p be an arbitrary point 
of Q = {(t,r) : r > 0,t*(r) < t < tnij)}. There is a unique outgoing radial 
null geodesic through p which by hypothesis extends back to the regular 
centre rather than the singularity. The inequality above tells us that any 
NRNG through p must cross the radial null geodesic from above, and so 
must precede this geodesic at all points q preceding p. Hence no NRNG's 
can extend back to the singularity • 

Now we consider the converse problem: Suppose that for a given mass 
function m(r) there are RNG's escaping from the singularity. Will there 
necessarily be NRNGs which also escape? To answer this problem we found 
it convenient to define new variables: 

x = i, y = R r f. 

Then a crucial point will be with how to control the derivative of x/y as 
r — > 0. We explain why this should be done before we actually prove our 
result in Proposition ahead. 
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We assume henceforth that there exist radial null geodesies emanating 
from the singularity. Choosing data at a point p with x\ p > 0, y\ p > 0, which 
corresponds to a future pointing outgoing null geodesic, we see that (x, y) 
lies on the positive right-hand branch of the hyperbola Ti l given by 



Clearly x/y > 1 along this branch. This says that the slope dt/dr in the r — t 
plane exceeds the slope of the unique future-pointing outgoing RNG through 
p. However, since p is an interior point of Q, we can find 5 sufficiently small 
such that the unique solution through pofx/y = l + 5 will also extend back 
to the singularity. (Notice that the proofs of existence of naked singularities 
for the radial case all involved showing that the slope dt/dr was smaller that 
the slope of some curve which we know to intersect the singularity.) Now 
on the hyperbola TCi, there is a continuum of choices of (x, y) which satisfy 
x/y < 1 + 5. Take any such pair and let this give the initial data for 
at some time So > 0. We follow the NRNG with these data as s decreases. 
Letting s — > + , we must have x/y > 1 for all s > 0, which says that 7 stays 
below the RNG through p. By our choice of data, 



to prove that 7 extends back to the singularity, it would suffice to show that 



This inequality says that 7 stays above the curve with x/y = 1+5, which 
we know extends back to the singularity. This would prove that the NRNG 
must also extend back to the singularity. The key to proving that we may 
proceed in this manner is to identify a subset of Q wherein the rate of change 
of x/y may be controlled. This is done by carefully analysing the relative 
magnitudes of the terms that arise in the equation governing this derivative as 
will be shown in the next proposition. (Note however that the demonstration 
of the existence of NRNG's emanating from the singularity does not quite 
follow the outline of the paragraph above: this paragraph is intended to 
motivate our attempt to control the derivative of x/y.) 




(38) 



{-){s )<l + 5; 



y 



(-)(*)< 1 + 5, 



s e (o,s ). 



y 
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Proposition 9 Let m{r) be such that there exist radial null geodesies ema- 
nating from the singularity. Let 

G(r,t)=R^. 

For r > 0, n > 0, 5 > define 

n[n,r ] = {(r,t) : < r < r , G(r,t) > * < t H (r)} 

and 

Q[5,n,r ] = {(r,t) : < r < r , yj - > 3^37(1 + S),t< t H (r)}. 

Then there exists r > 0, n > and 5 > such that for any p G Q[S, n, r ], an 
NRNG with initial data satisfying < 1+5 extends back to the singularity. 

Proof: Straightforward calculations show that (r,t) G Q[n, r ] if and only if 

t n (r) < t < t H {r) where 

3n ^ 2 y'rm 

n 3n — 1 c 3n — 1 rn r ' 

and that (r, £) G Sl[<5, n, r ] iff ts < t < t H where 

2 -1 
ts = t c - -cr m 

and a 1 / 3 = 2n(l + <5)/(3n - 1). 

The inequalities i*(r) < t n (r) < ts(r) < tuij-) may be shown to be true 
for sufficiently small r. This is easily done except in the final case, i.e. where 
m 1 (r) ~ — | / 9( i ^) 5 / 2 r 3 . In this case, the relevant inequalities will hold if n is 
chosen so that 

P > maxK = 3(3n - 1)(^)V 2 = 6 (^)(^) 3 } 

(In fact, the condition that (3 > n 2 arises from a later consideration, but we 
find it convenient to deal with it here). Rather fortuitously, we find that 
taking n\ = n 2 gives n\ = n 2 = and since we must have (5 > f3*, the 
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inequalities are verified. The inequality tg < t H requires that 5 > be 
chosen sufficiently small. 

Using the geodesic equations ( |3l3| - |3"5|) , we find that 



Now let p as defined be the initial point at time s = so of an NRNG and 
set the data such that (fp(so) < 1 + 5. By definitions, while the NRNG 
remains in Q[5, n, r ], we have ^(|) > 0. As s decreases, this geodesic cannot 
leave fl[5,n,r ] through t = tu (this would correspond to a future-pointing 
null geodesic exiting the trapped region across its past space-like boundary). 
Suppose that the geodesic reaches the boundary t = t${r). In the r — t plane, 
the geodesic has slope t' = |-R,r- But for s < s , while the geodesic is in the 
closure of Q[5, n, r ], we have (since x/y decreases as s decreases) 



Also, (1 + 8)R r \t=t 5 (r) < t's( r ) f° r sufficiently small r and 5. (As above, this 
is immediate except in the final case, where is necessitates the inequality 
(3 > n 2 .) Thus the geodesic cannot exit ft[S, n, r ] via this boundary, and 
therefore must extend back to the singularity • 

This then proves that if, for a given mass function m(r), RNGs escape 
from the singularity then NRNGs will also necessarily escape. The following 
theorem summarises the principal results of this section. 

Theorem 2 Given regular initial data for marginally bound spherical dust 
collapse subject to a condition which rules out shell- crossing singularities, 
there are non-radial null geodesies which emanate from the ensuing singu- 
larity if and only if there are radial null geodesies which emanate from the 
singularity. 

Finally, it is also of interest to note that there are always NRNG's which, 
starting from any point p with t < £#(?"), avoid the singularity in the past: 

Proposition 10 Let p lie in the region {(r, t) : r > 0,0 < t < %(r)}. If the 
initial data for an NRNG 7 through p satisfy (|)| p > §, then 7 cannot meet 
the singularity in the past. 




(39) 



t' < (1 + 6)R >r . 
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Proof: For < t < it is easily shown that 



G(r,t)<iy|. (40) 

Suppose that 7 meets the singularity. Then from above, it must do so in 
the region t* < t < and the condition that > x = t ^ — 00 in a 
neighbourhood of the singularity yields, using 

V < 2" 

But (69), (BO) and the initial condition at p yield t-(-)L < 0. Hence the 



condition xjy > 3/2 is maintained along the geodesic moving into the past, 
and so we obtain a contradiction. Hence 7 cannot meet the singularity in 
the past* 

Having proved the main results of this paper we will discuss some of their 
consequences in the next section. 



5 Conclusions and comments 

We have shown that the existence of radial null geodesies emanating from 
the shell-focussing singularity in marginally bound spherical dust collapse 
is a necessary and sufficient condition for the existence of non-radial null 
geodesies emanating from the singularity. Our interpretation of this result is 
that it expresses an aspect of stability of the naked singularity: the insertion 
of angular momentum into the system, at the level of the null geodesies, 
does not affect the outcome. This result can be interpreted as implying the 
robustness of the previous results, for radial null geodesies, with respect to 
the presence of angular momentum. One consequence of this is that all such 
naked singularities are equally 'malicious': previously, it has been noted that 
those which arise by the conditions (i) or (ii) of Theorem 1 being satisfied 
are gravitationally weak, while those which arise by condition (iii) of this 



theorem being satisfied are gravitationally strong [12]. However it is known 
that a non-radial null geodesic emanating from a central singularity always 
satisfies the strong singularity condition |Tj|, and so the distinction in terms 
of singularity strength vanishes when angular momentum is included. 

Our results on the relationship between the initial data and the causal 
nature of the singularity agree with previous results based on the study of 
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radial null geodesies. There are some differences: we have taken the initial 
data fi(r) to be C 3 , as opposed to analytic or C°°. The convenience of 
this assumption is that it is more general and allows a comprehensive and 
rigorous discussion of all possible cases. Can this level of differentiability be 
reduced? Some authors have studied this problem (in fact the more general 
problem where the collapse is not assumed to be marginally bound) using 
C 1 initial data (e.g. || §]; here a differentiability condition is imposed on 
the mass function m(r) which corresponds to /i(r) being C l ). However in 
these studies either a complete decomposition of the initial data space into 
regions giving rise to naked and censored singularities has not been obtained, 
or further restrictions on the differentiability are imposed. It seems to us that 
the critical role of the third derivative of /i(r) indicates that C 3 is as low as 
one can push the differentiability requirement while maintaining a complete 
description of the possible outcomes. 

As mentioned above, the approach used here to prove the existence of null 
geodesies emanating form the singularity is different from that of previous 
studies (but is more in line with the earlier studies @, 0)- One immediate 
advantage is that we can give independent verification of the previous results 
on radial null geodesies. We see other two distinct advantages to our (more 
qualitative) approach. Firstly, we do not need to assume the existence of a 
limit of R/r a (for some a > 0) in the limit as the singularity is approached. 
Assuming the existence of such a limit, or equivalently of y\{u)/u as u — > 
+ above, removes the need for the detailed analysis of Proposition 6: one 
can simply take the limit of both sides of (|32|) and use l'Hopital's rule to 
produce the quartic equation of Corollary 4. However, while it is almost 
certain that this limit exists for some of the null geodesies emanating from 
the singularity, it is by no means certain that this is the case for all such 
geodesies. Evidence that the limit need not exist arises from the contraction 
process in the proof of Proposition 6. If the limit did exist, one would expect 
this to be reflected in the bounding interval (k n , l n ) having zero length in 
the limit as n — > oo. In fact this does not happen, so while yi(u)/u is 
bounded within [lirn^oo k n , lim^oo l n ] in the limit as u — > + , there is no 
evidence that it possesses a finite limit. By way of illustration, we note that 
the approach relying on the existence of such a limit would not identify the 
(positive) solution 

y(x) — x(l + esin-) (0 < e < 1) 
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of the initial value problem 

y'(x) = ?L- -Jl- 4(1 --) 2 , 2/(0) = 0. 

Of course it may be that the regularity conditions imposed on the initial 
data only allow geodesies for which such a limit exists; however we feel that 
this is something which needs to be checked. Our fundamental point here 
is that the non-existence of this limit might not be a sufficient condition for 
the absence of a geodesic emanating from the singularity. 

A second advantage which we see is that the qualitative approach applied 
here lends itself very easily to the analysis of non-radial null geodesies. In fu- 
ture work we hope to show that this is also the case when one considers radial 
and non-radial time-like geodesies. This would then give a unified treatment 
of all causal geodesies in marginally bound spherical dust collapse. One clear 
disadvantage must be mentioned: We have dealt only with marginally bound, 
and not general, spherical dust collapse. We expect that the analysis used 
here can be carried over to the general case, but it must be acknowledged 
that this would not be done without some difficulty. In the general case, 
one no longer has the explicit solution (f4|) and consequently, identifying the 
important surfaces such as t = t*(r) becomes problematic. 

We comment as follows on the relationship between the initial data and 
the temporal nature of the singularity. One can argue that //(0) ^ 0, i.e. 
m'^O) 7^ is unphysical, as it represents a cusp in the initial density dis- 
tribution. One can also argue that all higher order odd derivatives should 
vanish at r = 0, i.e. that one should take /i(r) to be an even function of r 
|| The reasoning is that the polar co-ordinates (r, 6, 0) are singular at 
r = 0, and so one should move to a Cartesian co-ordinate system (x, y, z) in 
a neighbourhood of the origin. But then odd powers of r lead to non-smooth 
functions of (x,y,z). In particular, demanding C 2 behaviour in the Carte- 
sian co-ordinate system (it is difficult to see where higher derivatives would 
play a role in the physics, e.g. particle motion, of the space-time on and 
shortly after the initial slice) would rule out an r 3 term in /i(r) and enforce 
/x w (0) = 0. This leaves //"(0) as the only feature of regular initial data which 
has any input into determining whether the collapse results in a naked or 
censored singularity. Clearly, the generic case is /i"(0) < which leads to a 
naked singularity. 

Allowing the rather mild singularities generated by non-vanishing //(0) 
or yu'"(0), the picture becomes less clear. Obviously naked singularities will 
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arise as the generic case, but a more interesting question remains open: Do 
small perturbations of vacuum \x = and homogeneous \i = /i(0) initial 
data generically give rise to naked singularities? Consider, for example, the 
following situation. 

Fix the C 3 function mi(r), satisfying m'^O) = m"(0) = 0. Then we have, 
as in Section 3, m 1 (r) ~ — |/3(y-) 5//2 r 3 . Set initial data (fi ,m 1 (r)) and allow 
the collapse to proceed. Taking f3 < f3„, the singularity is censored. Now 
consider the initial data {p,Q,rhi{r) = mi(r)). Then 

Choosing p, sufficiently small, we can arrange that (3 exceeds the numerical 
value /3* and so the collapse will result in a naked singularity. The singularity 
will occur at time 

*o-g(y) >*o=g(y) ' 

Thus there is a tendency for 'smaller' initial data to result in naked rather 
than censored singularities. We note that this pattern reflects what is fre- 
quently found in studies of critical phenomenon: for example, in spherical col- 
lapse of a massless scalar field, Choptuik space-time, which contains a naked 
singularity, lies on the boundary between small/weak initial data which lead 
to dispersal, and large/strong initial data which lead to black hole formation 
jl9| . It would be useful to carefully analyse this issue as it arises in dust 
using a standard measure of initial data for general relativity p0| , |2l . 
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